Full perturbation solution for the flow in a rotating torus 
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We present a perturbation solution for a pressure-driven fluid flow in a rotating toroidal channel. 
The analysis shows the difference between the solutions of full and simplified equations studied 
earlier. The result is found to be reliable for low Reynolds number (Re) as was the case for a 
previously studied solution for high Re. The convergence conditions are defined for the whole range 
of governing parameters. The viscous flow exhibits some interesting features in flow pattern and 
hydrodynamic characteristics. 
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I. INTRODUCTION 

Fluid flow in a general curved pipe is known not to be 
one-dimensional since the pioneering work of Dean [l[ . A 
secondary flow develops when a fluid is driven by pressure 
gradient in a toroidal channel. In practice, this occurs in 
many applications: channels in industry, blood vessels in 
physiology, transport pipe systems, coolant pipe systems 
and others. In a considerable number of studies there has 
been a strong motivation for finding a solution for flows 
in a curved channel. 

A stationary solution of the Navier-Stokes equation 
in a curved pipe was first analyzed by the perturbation 
method pj. An analytical solution was obtained for a 
simplified equation valid at a low value of the Dean num- 
ber De = ReK 1 / 2 . This solution was compared with the 
results of direct numerical simulation for the wide range 
of governing parameters [2| . Later many authors consid- 
ered this problem using different approaches. Pressure 
drop, heat and mass transfers were studied in detail (for 
review see e.g. [2]). Non- Newtonian fluid flow in a curved 
pipe shows a specific behavior for different rheologic cases 
(e.g. Q and references within) . Magnetic field induction 
was experimentally studied in conductive fluid flow in 
toroidal channel j4j. 

In the present paper we derive a perturbation analyt- 
ical solution of the full equations taking into account all 
curvature effects. It is shown that the well known solu- 
tion [l| corresponds to the case Re 3> 1. In the general 
case the solution can not be parameterized by the Dean 
number only. The suggested solution properly describes 
the case of low Reynolds number (Re) and asymptotically 
approaches the known solution at high Re. For high vis- 
cosity of the fluid the inertial terms in the Navier-Stokes 
equation are relatively unimportant and the flow pat- 
tern is determined by a balance of viscous forces and the 
pressure gradients in the fluid. Such flows, called creep- 
ing flows, are of nearly the same importance in practice 
as the inertial flows (with Re ^ I). 



FIG. 1: The coordinate system in the toroidal channel. 



II. MATHEMATICAL MODEL 

A curved pipe is considered as a toroidal channel with 
the outer radius R c and radius of the inner circular sec- 
tion R. We use a coordinate system {p, ip, £}, where p 
and if are polar coordinates in the cross-section and £ is 
the linear coordinate along the channel (see Fig.[T|). The 
torus revolves around its main axis with a constant angu- 
lar velocity £1 (£1 > means co-rotation) . The flow of an 
incompressible fluid in the rotating channel is governed 
by the equations 



~dt 



v • Vv 



-2H x v 



Po 



i/Av, 



(1) 



V-v = 0, 



where v is the velocity, po is the fluid density and v is 
the kinematic viscosity. 

After introducing the vorticity lu = V x v pressure 
can be eliminated. The following Helmholtz equation is 
written as 



w + (vV) W - 



V)v = -2Vx(nxv)+i/Au. (2) 



Assuming that the flow is stationary and (^-independent 
Eqs. ©-(O can be reduced to the equations for two 
scalar functions: stream velocity and stream function 
ip which is defined in such way that 
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The no-slip boundary conditions 



dp 



(4) 



together with the regularity condition at the centerline 
must be satisfied. Following arguments in [f| we assume 
the pressure to be a linear function of £. Therefore the 
pressure gradient is constant vector along channel. 

Choosing V as the characteristic velocity, the inner 
radius R as the characteristic length and Q as the char- 
acteristic pressure gradient we introduce dimensionlcss 
variables p* = p/R, (* = (/R, v* c = v c /V, ip* = 
ip/(VR), p* = pR/K. The following governing di- 
mensionless parameters will be used: curvature ratio 
k = R/R c , Reynolds number Re = (VR)/v, dimension- 
less angular velocity F = (flr c )/V and pressure parame- 
ter G — (QR 2 ) / (pqvV) which is a ratio of the viscous and 
pressure forces, r — — p cos ip) is the dimensionless 
distance from the torus axis. 

Thus, the following set of equations describes the prob- 
lem (the superscript '*' is omitted everywhere below) 
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p dip 

while the derivation of the curl in these coordinates leads 
to the expression for stream function: 



v = D v ip 
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I dip 



v v = -D r ip = 
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Other differential operators are 
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The effect of curvature can be obtained after significant 
simplification of the governing equations {5} and ^ . The 
idea of Dean [l[ was to renormalize the stream function 
by a factor of Re and combine factors Re k 1 / 2 into the 
so called Dean number. The approximation k <C 1 is ap- 
plied after that. Terms due to curvature remain only in 
the first two advective terms and the Coriolis force term 
in Eq. ([6]). Limit k <C 1 and fixed value De corresponds 
to limit Re 3> 1. A perturbation solution of these simpli- 
fied equations was reconsidered and compared with the 



results of direct numerical simulations in Q. We solve 
the equations ([5]) and © in general. This means that 
the viscous term corresponding to the toroidal Laplacian 
operator C (fTUj) is not replaced by a cylindrical one ©. 
The term corresponding to the Coriolis force in Eq. ([5]) 
remains. All curvature contributions are kept in the ad- 
vective terms and pressure gradient. 



III. PERTURBATION SOLUTION 

Solution of Eqs. (O and (JSJ) can be expanded in power 
series in a small parameter k 

V C = U (0) (P, <P) + KW (1) (/9, ip) + K 2 V^(p, tp) + (11) 

ip = V> (0) (p, ip) + ™/> (1) (p, <p) + t V 2) (p, <p) + ■ ■ • (12) 

Expanding all operators ([7|)- (fT0|) in the same series we 
can write the zero-order approximation of Eqs. (JSJ) and 
© as 

f,M ,/.(o)\ _ G_ , _L A „(o) 
V J-Re + Re Ac ^ ' 
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(13) 



where J [A, B) is the Jacobian of the functions A and B 
with respect to p and ip: 

1 fdAdB dBdA s 



This gives a solution corresponding to cylindrical geom- 
etry, namely to the Poiseuille flow 



,(°) = ^(l-p 2 ),^)=0. 
The first-order approximation leads to the system 

J 



(15) 
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(17) 



and its solution gives the first-order correction 
3G cos ip Ia 2 \ 
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p(l- P 2 )(p 4 -3p 2 + 3), (18) 
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This solution is different from the solution obtained by 
(we denote the latter by superscript 'c'). The first-order 
solution has an additional term 



i(e) _ K 3Gcos^ p (1 _ p2)+Q (k2) (20) 



U C - 16 

has no difference from that in 0]. The corrections 
of the second-order approximation and can be 
derived in a similar way. For the stream function we find 

oP 2 (l - P 2 ) 2 sm2w , 

ib - i/; (c) = — k — — 150GReF+ 

y Y 92160 v 

+G 2 Rc(56- Up 2 )) +0(k 3 ) . (21) 

We note that K and Re cannot be combined (for instance, 
in the Dean number) for parametrization of the present 
perturbation solution, both should be treated as inde- 
pendent parameters. 

IV. CONVERGENCE 

It is clear that the perturbation solution is valid only 
at specific values of the governing parameters Re, F, G 
and k and the series must converge. The conditions on 
the governing parameters have been defined in order to 
satisfy the uniform convergence of the maximum value 
for each order of approximation. 

All terms in the expressions for «W and ip® have the 
form (GRe) J (FRe)' c /Re. So, for convergence analysis we 
use for convenience the following values: G' = GRe and 
F' = FKe which are the only factors that influence the 
convergence. The maximal powers of G' and F' included 
in t)W grow as 2i + 1 and i accordingly. However, these 
powers can be reached only in different terms such as 
q/2i+i anc j qii+i p'i /R e . So the radius of convergence 
should be 0(G'~ 2 +F'~ 1 ). We estimate the radius of con- 
vergence k c with the following technique. A supposition 
of uniform convergence neccessarily implies a monotonic 
decrease of «W majorants in (flTj) . In order to know this 
we evaluate terms up to the 5th order, find the global 
maxima in the circle p < 1, and calculate largest k, for 
which the maxima do not increase. 

Thus we obtain function of G' and F' (see 

Fig. [2]), rather complex but holding some nice asymp- 
totics. For large F' we have obtained a dependance of 
the following kind: a + b \F' — c|. As pattern fitting has 
shown, the coefficients a, b and c grow linearly in G' . The 
slope of the infinite branches grows as b ~ G' . 



V. RESULTS 

Detailed analysis of the solution was done in the earlier 
studies. Where the attention was centered on the case of 




FIG. 2: The dependence of convergence radius on F' for G' 
5. .705 (top) and for G' = 5. .21000 (bottom). 




FIG. 3: « = 0.5, Re — 1, F — 0. Comparison of solutions. 

(a) v c , (b) vf\ (c) V, (d) V (c) - 



large Re. Our new finding deals with the case of small Re 
when the viscous term is comparable or larger than the 
convective one. This difference arises when new terms in 
the solution (|20H2ip dominate, i.e. when Re is small. 

In the case of F = the result is shown in Fig. [3] One 
can see that the shift 5 m of the stream velocity maxi- 
mum due to a centrifugal force directed outward from 
the torus axis is not valid at small Reynolds number. In 
fact for k > 0.3 there is a shift only in the opposite direc- 
tion (inwards) because of viscous stresses in the region 
of convergence. This shift is mainly due to the difference 
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FIG. 4: Shift S m of the maximum of stream velocity vs Re at 
k — 0.1. Dashed line stands corresponds to the solution v^ c \ 




FIG. 5: k = 0.5, Re = 1, 
tions. (a) ip, (b) ip (c) . 



F = —0.63. Comparison of solu- 



(|2"U)) . It has a dependence on Re shown in Fig. [U 

When the torus rotates (F ^ 0) the Coriolis force gives 
rise to additional vortices in the cross-section. In the 
counter-rotating case such a vortex can act against the 
centrifugal vortex. This produces a four-vortex picture 
(see Fig. [5]) . Again we see that the flow pattern is differ- 
ent for a low Reynolds number. The vortex correspond- 
ing to the Coriolis force arises at the boundary while tp^ 
starts to grow in the center. 



Axial shear-stress along the channel r = — 



dif- 



P =i 



fers from that presented in [2J. The first order residual 
is 3/16Gk (3p 2 — l) cost^. It does not change the pres- 
sure drop but produces a strong variation (about 40%) 



of friction at the boundary (see Fig. [6l thick curves). 
There is a similar difference in the azimuthal shear-stress 



dp 



The second order residual is 



GRe(13G-25)K 2 sin2ip 
3840 ' 



(22) 



This work shows that solution of the full governing 
equations reveals some specific features for flow in a 
toroidal channel: a shift of the maximum of the stream 
velocity toward the inner axis, the appearance of a second 
pair of vortices at the internal boundary, an additional 
</?-dependance of stresses. These are well pronounced at 
low values of the Reynolds number. The solution for a 
giving curvature k asymptotically approaches the known 
solution Q at high Re. 




FIG. 6: Axial shear-stress r (bold lines, tick labels at left) and 
azimuthal shear-stress r v (thin lines, tick labels at right) de- 
pendencies on angle: Dashed lines corresponds to the solution 
in 01. k= 0.4, Re = 10, F = -0.49. 
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